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Abstract 



. By using the thermodynamic theory of irreversible processes and Ein- 

stein general relativity, a cosmological model is proposed where the early 
universe is considered as a mixture of a scalar field with a matter field. 

■ The scalar field refers to the inflaton while the matter field to the classical 

particles. The irreversibility is related to a particle production process at 
the expense of the gravitational energy and of the inflaton energy. The 
particle production process is represented by a non-equilibrium pressure 
in the energy-momentum tensor. The non-equilibrium pressure is prepar- 
es tional to the Hubble parameter and its proportionality factor is identified 

with the coefficient of bulk viscosity. The dynamic equations of the in- 
, flaton and the Einstein field equations determine the time evolution of 

the cosmic scale factor, the Hubble parameter, the acceleration and of 
the energy densities of the inflaton and matter. Among other results it is 
shown that in some regimes the acceleration is positive which simulates 

■ an inflation. Moreover, the acceleration decreases and tends to zero in the 

instant of time where the energy density of matter attains its maximum 
value. 

PACS: 98.80.-k; 98.80.Cq 



1 Introduction 

Cosmological models are among the most important formulations that can be 
put under analysis when we combine the thermodynamics of irreversible pro- 
cesses with Einstein's general relativity pi |l H, m H, @, m, @, 0, |uj|. Using 
these theories the different eras of our universe can be modeled; the subject 
received great attention: the starting point was to consider the hypothesis of 
homogeneity and isotropy in the form of the Robertson- Walker (RW) metric. 
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More specifically, the early universe physics is of great interest; the so-called 
inflationary formulations focus on that era, when - according to these models - 
a rapid expansion of space-time would occur. Using these ideas it is possible to 
construct a non-equilibrium scenario, where a non-equilibrium pressure term is 
included in the energy-momentum tensor; this introduces the possibility of cre- 
ation of particles in the early universe by means of the thermodynamic theory of 
irreversible processes. In the 3+1 case the inclusion of this term was proposed 
by Murphy || who found a solution that corresponds to a simple expansion. In 
the inflation driven models a fundamental field, the inflaton, simulates the early 
universe with an equation of state that will adjust the form of the potential 
and will affect the temporal evolution of fundamental quantities like the cosmic 
scale factor and the universe's energy density. 

In this work we start with the Einstein field equations, the inflaton dynam- 
ics (in a curved space-time), the conservation law for the energy-momentum 
tensor, the barotropic equations of state of the inflaton and matter fields. The 
non-equilibrium pressure is considered as a constitutive quantity - within the 
framework of ordinary (also known as first order Eckart) thermodynamic theory 
- which is proportional to the Hubble parameter and whose proportionality fac- 
tor is identified with the coefficient of bulk viscosity. From these equations we 
determine the time evolution of the cosmic scale factor, the Hubble parameter, 
the acceleration and of the energy densities of the inflaton and matter. The evo- 
lution equations for these quantities are function of three parameters which are 
identified with the cosmological constant and two coefficients that appear in the 
barotropic equations of state for the pressures of the inflaton and of the matter. 
It is shown that among several possible regimes the acceleration is found to be 
a positive quantity which simulates an inflation. Moreover, in those regimes the 
acceleration decreases and tends to zero in the instant of time where the energy 
density of matter attains its maximum value. 

The manuscript is structured as follows. In Section II we introduce the bal- 
ance equations for the particle four-flow, energy-momentum tensor and entropy 
four-flow and use the Gibbs equation to identify the non-equilibrium pressure 
with the process of particle production. We discuss the inflaton equations in 
Section III and obtain the evolution equation of the energy density of the in- 
flaton as a function of the cosmic scale factor. In Section IV we introduce the 
Einstein field equations and find a second-order differential equation for the evo- 
lution of the cosmic scale factor which is a function of three above mentioned 
parameters. In Section V we search for the evolution equations of the cosmic 
scale factor, Hubble parameter, acceleration and of the density energies of the 
inflaton and matter for given values of the three parameters and initial condi- 
tions for the cosmic scale factor and Hubble parameter. Finally, in Section VI 
we consider the non-equilibrium pressure as a variable within the framework 
of extended (also known as causal or second-order) thermodynamic theory and 
discuss the behavior of the solutions. Throughout this paper units have been 
chosen so that the speed of light c and Planck constant fi are equal to one. 
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2 The Balance Equations 



According to the cosmological principle the universe is spatially homogeneous 
and isotropic so that it looks the same to all observers. These assumptions imply 
that the universe can be described by the Robertson- Walker metric whose line 
element in a space-time, characterized by the metric tensor with signature 
(+ — ), is given by 



ds z = dt 1 - a(ty 



dr 2 



1 — kr 2 



+ sin 2 6dip 2 



(1) 



In the above equation a(t) is the cosmic scale factor which is an unknown func- 
tion of the time t and k may assume the values 1, 0, -1 that correspond to 
closed, flat and open universes, respectively. 

We consider that the universe can be modeled as a mixture of two con- 
stituents, namely a scalar field and a matter field, in interaction via the gravita- 
tional field represented by the cosmic scale factor a(t). There are other ways of 
modeling this interaction and for more details one is referred to Yokoyama and 
co-workers , Bercra , Billyard and Coley a nd the references therein. 
The scalar field refers to the so-called inflaton Jl4|, |15[ and the matter field de- 
scribes the particles classically. The mixture is characterized by the fields of 
particle four- flow 7V M , energy- momentum tensor and entropy four-flow 
whose balance equations read 

A^ ;M = n£, T^ v . v = 0, S% > 0, (2) 

where n denotes the particle number density while £ is the particle production 
rate. 

In a homogeneous and isotropic universe the decompositions of the particle 
four-flow, energy-momentum tensor and entropy four-flow with respect to the 
four- velocity (with U^U^ = 1) read 

N" = nU l ", r"" = (p+p + w)U"U u - (p + wW, S fl = nsU fl , (3) 

since in this case the stress tensor, the heat flux and the entropy flux vanish. 
In (^) p, p and s represent the energy density, the pressure and the entropy 
per particle of the mixture, respectively. The term w denotes the dynamic 
pressure which refers to a non-equilibrium pressure; this quantity also represents 
the production of particles phenomenologically p6[ . For this particular case 
the Eckart and the Landau and Lifshitz decompositions (see for example p7|) 
coincide. 

Insertion of the representations (|^) into the balance equations iV M :/J = nS, 
UpT» v .y = 0, and 5% > lead to 

n + n6 = nS, p + (p + p + zu)Q = 0, ns + nsS > 0, (4) 

which represent the balance equations for the particle number density, energy 
density and entropy density, respectively. In these equations we have introduced 
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the usual notations: 9 = U^-^ for the expansion rate and h = U^n-^ for the 
covariant derivative along U^. 

In terms of the the covariant derivative along the Gibbs equation reads 

ns=±(p-P±Pny (5) 

where T is the absolute temperature. If we eliminate from the Gibbs equation 
the derivatives of the particle number density h and of the energy density p by 
using the two first balance equations given in (|^) we get 

ns + nsE = — (— vj<d — npJS) > 0, (6) 

where p, = (p + p)/n — Ts is the chemical potential. 

An adiabatic process is characterized by the requirement that s — so 
that we can infer from (^|) the inequality nsS > 0. As was pointed out by 
Prigogine et al the process of particle production leads to a non-symmetric 
relationship between space-time and matter since it refers to an irreversible 
process of particle production due to the gravitational energy. Moreover, we get 
from (^) when s = that 

w = -{p + p)-. (7) 

Hence, the dynamic pressure w can be identified with the particle production 
rate E. This result was first obtained, to the best of our knowledge, by Zim- 
dahl |§. 

We write the energy density and the pressure of the mixture as a sum of two 
terms that describe the scalar and the matter fields, e. g., 

p = p<t, + p m , P = P4>+p m , (8) 

where p$, p$, and p m , p m represent the energy density and the pressure of the 
inflaton and of the particles, respectively. In this case the energy-momentum 
tensor for the mixture reads 

= fa + Pm+Pt + Pm + m^U" - (pt +p m + m)g» v . (9) 

From this point on we shall consider a comoving frame where the four- 
velocity is given by (U^) — (1,0). In this circumstance the expansion rate is 
given by O = 3H - where H — a/ a is the Hubble parameter - and the balance 
equation U^T^.^ = can be written as 

+ Pm + 3H(p0 + Pm + P4> + Pm + w) = Q. (10) 

In a comoving frame the dot reduces to a differentiation with respect to the 
time. 
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3 The Inflaton Equations 

Modern inflationary theories require the presence of the inflaton. This hypo- 
thetical particle is represented classically by a scalar field of the same name and 
the motivation behind these ideas can be found in the works JTJ, |f5[ [fj| |T| . 
The corresponding Lagrangian density (in a generic curved space-time) for the 
scalar field (j)(x^) is written as 

£ = ^<OT-V(», (ff) 

where V (</>) is the potential density of the field. From the Euler-Lagrange equa- 
tions of motion one can obtain the dynamics of a homogeneous scalar field 

+ 3^+ ^M = 0. (12) 

a<p 

If we consider the inflaton as a perfect fluid with an energy-momentum 
tensor of the form (T£ v ) = diag(p0,p^,p0,p^,), we can equate this formula to 
the expression coming from Noether theorem 

T£" = d»<t>d v <t>- CgT, (13) 

and obtain the equations for the energy density and for the pressure in terms 
of the homogeneous scalar field, which are given by 

P4> = l^ + V(4), P4> = \p-V{<t>Y (14) 

We differentiate (jTJ)i with respect to the time, make use of ( fl2|) and get the 
following evolution equation for the energy density of the inflaton 

Pt + SHipt+p^) =0. (15) 

Equations ( |Io| ) and ( |l5|) show that the energy densities of the inflaton and 
of the matter decouple so that we can write 

p m + 3H(p m + p m + w) = 0. (16) 

Now we follow the work of Ratra and Peebles |1| and suppose that the 
pressure of the inflaton is connected with its energy density according to the 
barotropic equation of state 

Pcf, = (v — l)p0, where < v < 1. (17) 

As we said before, in this model the homogeneous scalar field does not interact 
with other non-gravitational fields |19| . 



We insert (17) into (|15[) and get by integration a relationship between the 



energy density and the cosmic scale factor 

(18) 
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where p° and ao are the values of the energy density of the inflaton and of the 
cosmic scale factor at t = (by adjusting clocks), respectively. The relationship 
between the potential density of the scalar field and the cosmic scale factor can 
be obtained from Q, (jl7|) and (|l|) yielding 

The above result is in agreement with the work of Ratra and Peebles jl9| as 
expected. 

Once the time evolution of the cosmic scale factor is known it is possible 
to obtain from (|l8|) and ([lj]) the time evolution of the energy density and of 
the potential density of the scalar field. The determination of the cosmic scale 
factor from the Einstein field equations will be the subject of the next section. 



4 Einstein Field Equations 

From Einstein field equations 

Rfjtv — ^Rgiiv = ~-^T IJjU , (20) 

where mp = 1/vG denotes the Planck mass, R^ v the Ricci tensor and R the 
curvature scalar, one can get (working with the RW metric) a system of equa- 
tions that reads 



47T 

3^1"' 



+ p m + 3p0 + 3p m + 3ra7), (21) 



k 8tt 



H -\ J = o 2^P<t> + Pm)- (22) 



p 

The system of equations (|l0|), ( pl| ) and ( |2^ ) is not linearly independent, since 
the differentiation of (22) with respect to the time, taking into account the 
conservation law (|Io| ) leads to (|2l|). 

Let us first analyze the case that correspond to a false vacuum where v = 
so that prj, = —p<f> = —p\- In this case we have only the scalar field, i. e., the 
matter field is absent (p m = 0) and it follows from (P2|) 

- = -^—t = X, (23) 

where x can De identified with the cosmological constant A = 3%. Equation 
( p3[ ) leads to the de Sitter exponential solution a = ao exp (^/%t) if one neglects 
the term k/a 2 during the rapid expansion. 

For v ^ we can get from (^2|) together with ( |Io|) 

fc 47T 

-H" = ^ - ^2"(P0 + Pm +P0 +Pm + ^)- (24) 



G 



In order to determine the time evolution of the cosmic scale factor from 
24]) one has to know the relationship between the pressure p m and the energy 
density p m of the matter and the constitutive equation for the dynamic pressure 
w, since p$ = {y — l)p^ and the energy density of the scalar field p^ is given by 
Here we follow the works and write 



Pm = (7 - 



1<7<2, 
rj = a(/90 + p m ). 



(25) 
(26) 



Equation ( |25| ) is the barotropic equation of state for the matter and some values 
for 7 are: a) dust 7 = 1; b) radiation 7 = 4/3; c) non-relativistic matter 7 = 5/3 
and d) stiff matter (or Zel'dovich fluid Q) 7 = 2. Equation (|26| ) relates the 
dynamic pressure w with the Hubble parameter H . The proportionality factor 
77 is the coefficient of bulk viscosity which is supposed to be proportional to the 
energy density of the mixture p = p^ + p m . Moreover, a is a constant. 

The evolution equation for the cosmic scale factor is obtained from ( pi| ) 
together with (H) and ([||) yielding 



H 



{v - i)x ( 



a J 



3u 



(7- 





Ml 


(h 2 h 




a 2 ) 



(27) 



The above equation is a function of four parameters - namely v, \i 7 an d ot - 
for a given value of k. We can express one of these parameters as a function of 
the others. Indeed if we consider that at t = (by adjusting clocks) we have 
H(0) = H , p m (0) = and a(0) = a it follows from (||) and @ 



tt2 ^ 



9H oX a 2 

Now we make use of the relationships ( p8| ) and get from ( 27 ) 

Zv X - 2fc 



(28) 



H = 



a 2 (\- k) 



3X 



3, _ s X 
2 X — k \ (i 



H H 



a 2 (x - k ) 



(29) 



by introducing the dimensionless quantities: a) a new Hubble parameter H = 
H/Hq; b) a new proper time t = tHo; c) a new cosmological constant x = X a o 
and d) a new cosmic scale factor a = a/a^. Equation (|2^) is a second-order 
differential equation for the cosmic scale factor which is a function of three 
parameters: the cosmological constant Xi an d the two coefficients that appear 
in the barotropic equations of state for the inflaton v and matter 7. 
In terms of the dimensionless quantities the energy densities read 



3w 



x-k 



a 2 (x~- k) 



3// 



(30) 
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The above equations connect the evolution of the energy densities with the 
evolution of the cosmic scale factor. In the next section we analyze the time 
evolution of these quantities for flat, closed and open universes (k — 0, ±1) and 
for given values of the parameters is, x and 7. 

5 Results and Discussions 

In order to obtain the time evolution of the cosmic scale factor and of the energy 
densities from equations ( ^9|) and (^) one has to specify two initial conditions, 
since (|2^) is a second-order differential equation for a(t). Here we choose that 
at the instant of time t = (by adjusting clocks) the energy density of the 
inflaton has its maximum value while the energy density of the matter attains 
its minimum value, namely p$(0) — p\ = 1 and p m (0) — 0. According to 
( |30| ) these initial conditions are equivalent to fix the value a(0) = 1 for the 
cosmic scale factor and H(0) — 1 for the Hubble parameter. There is still 
much freedom to find the solutions, since (|2^) and (|3(]) do depend on the three 
parameters v, x an d 7 for each kind of universe, i. e., closed, open and flat. 
Among several possible regimes which can be found by choosing different values 
for the parameters v, x and 7 we fix our attention to those values which simulate 
an inflation - i. e., where the acceleration is a positive quantity - and where the 
energy of the matter attains its maximum value at the instant of time where the 
period of the inflation ends - which corresponds to a vanishing acceleration. For 
a closed universe (k = 1) one set of values that satisfies the above mentioned 
conditions are v = 0.9, 7 = 1.9 and x — 2.0. In figures 1 and 2 we have also 
chosen the same values for a flat (k = 0) and open (k = —1) universes in order 
to compare the different solutions. 

In Fig. 1 it is plotted the cosmic scale factor a, its velocity d and acceleration 
a as a function of time for closed, flat and open universes. We infer from the 
curves that an open universe evolves more rapidly than the two other types of 
universe, while a closed universe evolves more slowly. The same conclusion can 
be draw out for the acceleration, since the curves show that the inflation period 
of an open universe is the largest following the flat and the closed universes. We 
have plotted in Fig. 2 the energy densities of the inflaton and of the matter as 
function of time. We note from this figure that the energy density of the matter 
grows more rapidly with respect to the time for an open universe following 
the cases of flat and closed universes. Moreover, we can infer that the energy 
density of the inflaton decays more slowly with respect to the time for a closed 
universe following the cases of flat and open universes. Once the maximum of 
the energy density of the matter is attained other models should be considered 
for the universe, since the universe reheats and follows the standard big-bang 
model. 

One may wonder from Fig. 2 why the sum of the energy densities of the infla- 
ton and of the matter is not a constant. This can be explained as follows: in the 
presence of gravitational fields the energy-momentum tensor of the matter alone 
does not lead to a conservation law (see for example Landau and Lifshitz p2[ 
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and Dirac p3[), since one has to include the energy-momentum pseudo tensor 
of the gravitational field Tq 1 which is given by 



9 \ l a/3 1 <jt 1 aa 1 Bp)9 



, (31) 



where rg, CT are Christoffel symbols and g = — dct(<7 AI „). 
For a RW metric with k = we get 



TlOO _ 

1 G — 



Snip ( a 



8tt \a 



(32) 



hence giving the "expected" relation T°° + T°° + = thanks to T°° = 

p m and d22) . One may say that the energy density of the matter increases 



at the expense of the inflaton and gravitational energies. This result must 
be handled with care since the energy density of the gravitational field is an 
intrinsically non-covariant quantity p3, E3[ . 



6 Final Remarks 

In this work we have analyzed the irreversible processes within the framework 
of ordinary thermodynamics instead of extended thermodynamics (see for ex- 
ample [| [?], |8|, |^]). The reason is that - apart from causality - there exist 
some thermodynamic systems in which ordinary and extended thermodynamics 
lead to the same results, like the problems concerning the propagation of low 
frequency waves in fluids and those related to the flow and the heat transfer 
in fluids in the continuum regime (see for example |2^| ). When the behavior 
of the thermodynamic fields are not smooth enough, there exist differences be- 
tween the solutions of ordinary and extended thermodynamics. Here we confirm 
that both descriptions can lead to the same behavior of the solutions and for 
this purpose let us proceed to analyze the same problem by using extended 
thermodynamics . 

In extended thermodynamics the dynamic pressure w is no more given by 
the constitutive equation (p6| ) but it is a variable whose evolution equation in a 
linearized theory reads (see for example |l7], |24| ) 

zu + tvj = —3rjH, (33) 

where r is a characteristic time. The evolution equations for the heat flux and 
pressure tensor are not consider here, since we are dealing with a homogeneous 
and isotropic universe. Hence we have to solve the system of evolution equations 



for the cosmic scale factor (|24|) and for the dynamic pressure (33). 

According to the kinetic theory of relativistic gases (see for example p"7[) 
the characteristic time which appears in the evolution equation for the dynamic 
pressure is given by t = i]/p, where p is the pressure of the mixture. Hence 
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the system of equations (£ 
quantities as 



and (J33|) can be written in terms of dimensionless 



H- 



a(x - k)w + 



(v-l)x ( 1 
X-k \a 

w + 3aH ( H 



3// 



1 [ H 



k 



a?{x - k) 



0, (34) 



a 2 (x- k) 



a 



3w 



+ ( 7 -l)( X -fe) iJ 2 + 



a 2 {x~ k) 



= 0, 



(35) 

to/(3w^ 2 ), 



by introducing the dimensionless quantities a = ai/o and 
apart from those defined in Sec. IV. 

As we have pointed out in previous sections there exist several possible 
regimes that can be found by solving the system of equations ([m|) and (|35|), 
since it is a function of four parameters v, Xi 1 an d ex. If one chooses the 
following initial conditions a(0) = 1, H(Q) = 1 and wifi) = one may obtain 
for (say) v — 0.9, x — 2-0, 7 = 1.9, a = 0.1 and k = 1 a solution for the fields 
in extended (causal) thermodynamics which have the same behavior as those 
represented in Figs. 1 and 2 by using ordinary (Eckart) thermodynamics. 
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Figure 1: Cosmic scale factor a, velocity a and acceleration a vs time t for 
open k = — 1 (dot-line), flat k = (dash- line) and closed k = 1 (straight-line) 
universes. 
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Figure 2: Inflaton energy density and matter energy density p m vs time t for 
open k = — 1 (dot-line), flat fc = (dash- line) and closed k = 1 (straight-line) 
universes. 
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